Thermodynamics of compact-star matter within an Ising approach 
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7— i ■ In the formation and evolution of compact stars, nuclear matter explores high thermal 

^ . excursions and is the site of intense neutrino emission. Neutrino transport as well as struc- 
tural properties of this matter depend on the presence of inhomogeneous phases (named 
"pasta" phases), which are the result of Coulomb frustration of the Liquid-Gas phase 

y—{ ; transition. We take into account charge fluctuations by employing a frustrated lattice-gas 
model to which we impose a neutrality constraint by the addition of an homogeneous 
background of charge, representing delocalised electrons. Within this schematic model 
we highlight a generaic feature of the phase-transition phenomenology: the temperature 

i— h \ interval where pasta phases are formed is enhanced by Coulomb-frustration effects. This 
result is at variance with the behaviour of frustrated ferromagnetic systems as well as 
hot nuclei and mean-field approaches. Moreover, the region of phase coexistence is not 
J> ■ found to end upon a critical point, indicating that no critical opalescence can occur in 

^ . compact-star matter. 

£ : 

1. Introduction: stars and atomic nuclei 

Nuclear matter present in compact and hot astrophysical objects generated in a gravita- 
tional collapse suffers from high thermal excursions. In the accessible range of temperature 
and densities, nuclear matter presents phase transitions characterized by density fluctu- 
ations. Because of the presence of the long range Coulomb interaction, those charge 
fluctuations give rise to so called "frustration effects" [H]. Under these conditions, the 
system organises in inhomogeneous phases, named "pasta phases" due to their unusual 
topologies [ 2,. The evolution of these structures with density is expected to be con- 
nected to the properties of the inner neutron-star crust. The survival of pasta-phases at 
finite temperature has been inferred from recent molecular-dynamics simulations [Ej. As 
these complex structures are related to coherent neutrino scattering their evolution 
with density and temperature or, more generally, the underlying equation of state, be- 
comes a key quest for exploring several interconnected astrophysical phenomena, like the 
dynamics of supernova explosion and the thermodynamics of proto-neutron-stars cooling. 

A tempting connection can be searched with the physics of atomic nuclei, which are 
experimentally accessible. Comparable conditions of temperature and density can be 
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probed in rather violent processes like nuclear multifragmentation. However, compact- 
star matter differs from normal nuclear matter due to the presence of electrons. They 
constitute an incompressible degenerated gas, which establishes a condition of global 
charge neutrality at large scale. 

To understand the phenomenology of phase transitions in compact stars, we consider 
an Ising model. Indeed widely used mean-field approximations [ El El El are known to 
fail to describe critical phenomena in 3D. The Ising model is specifically suited for the 
analysis of charge fluctuations by exact calculations, since phase transitions are universal 
phenomena. Moreover, it was largely employed in the study of ferromagnetic systems 
subjected to frustration as well as for nuclear matter where, in the form of the 
lattice-gas model, magnetic spins are replaced by site occupations. 

2. Ising analogue to compact-star matter 

We construct a cubic periodic lattice composed of V sites, each characterised by a po- 
sition i and an occupation number = 1 or 0, to indicate the presence or the absence 
of positive charge, respectively. This distribution of positive charge simulates the distri- 
bution of nuclear-matter charge (neutron and protons without distinction are considered 
as nucleons with an effective charge Z/A). The strongly incompressible gas of electrons 
is represented by a uniform distribution of negative charge [ imposing a strict con- 
dition of neutrality. In any site we find therefore an effective charge q\ = n\ — n, where 



n 



J2\ n slV is the charge per site of the background of negative charge. 



The schematic Hamiltonian H n+ c = Ti n + Tic with 
W» = f£Vj, ^c = yEf^ (1) 

is introduced to study the interplay of nuclear-like (K n ) and Coulomb-like forces (Tic)- 
X/ij is a sum extended over closest neighbors, and is the distance between sites i and j. 
The short-range and long-range interactions are characterized by the coupling constants 
e and Ae = dhep^x 2 respectively, where po is the nuclear saturation density, and x is 
the proton fraction. The ratio of the two coupling constants A measures the strength of 
frustration. In the sequel, x is fixed to 1/3, a value expected for proto-neutron stars. Tic 
can be rewritten in terms of occupation products times the geometric constants Cy as 

?*c= yEn^Cy. (2) 

ij 

Translational invariance, which is ensured by imposing the distance rij to be the shortest 
between i and j in the periodic space, imposes Cy = 0. Hence, Cy can be calculated 
by adopting periodic boundary conditions in the form 

1 v 

Cij = Aj - T>, with V = -Y i D 0V , (3) 

i' 

where = if i = j, and = | r i j | 1 otherwise; T> is calculated with respect to any 
position 0. 
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Tic preserves all the characteristic symmetries of the Ising model and the Coulomb 
energy is therefore invariant for lattice-occupation reversal. We should also point out 
that, despite some tempting similarities, the correspondent of eq.Q in spin variable 
Si = rii — 1/2 and p = M./V + 1/2 with the magnetization M. = Ya «i is not identical to 
the Hamiltonian of the frustrated Ising ferromagnet [EI], defined by 

Wi* = f f>j + y X>jAi- (4) 

u u 

Indeed, Ti-p e and Ti n+ c are related by 

H Fe = H n+C + -^VM 2 + 3eM - jV, (5) 

This connection includes a term in A4 2 , scaling with V 2 . Such term deeply modifies the 
thermodynamics in a finite system, and imposes M. = at the thermodynamic limit 
for Ti.F e . The strict constraint M. = on the order parameter M. leads to a substan- 
tial modification of the thermodynamics of the system[EI|. The phase diagram of the 
frustrated ferromagnet [E3> where , as in nuclei, the Coulomb interaction reduces the 
phase-coexistence region (i.e. decreases the limiting temperature), cannot be directly 
exported to our system 7i n+ c- 

3. Simulation in the multi-grandcanonical ensemble 

For the numerical simulation of the H n +c system, in order to take into account the long- 
range Coulomb interactions in a rapidly converging form, the finite lattice is repeated in 




Figure 1. Example of a Metropolis calculation (/i = p c ). Density distribution for a 
fixed value of /3c and a series of temperatures, each one associated to one spectrum. 
The ensemble of explored temperatures and densities defines the phase diagram. The 
ridge of the distribution defines three branches joining at the limiting temperature. The 
phase-coexistence region is delimited by the gas (low density) and liquid (high density) 
branches. 
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all three directions of space a large number R of times. Each site i has R replicas of itself, 
each one displaced from i of a vector mL, where m has integer components and L = \/V. 

To treat a frustrated system, described by the interplay of several interactions, the 
multi-(grand) canonical ensemblef E] is specifically suited. It consists in associating a 
Lagrange multiplier to each component of the Hamiltonian that we want to treat as a 
separate observable. In particular, it makes possible to construct one single phase diagram 
for both neutral and charged matter by associating one Lagrange multiplier, (3 n , to the 
nuclear energy (related to TC a ) and another, /3 C , to the Coulomb energy E c (related to 
T~Cc)- The multi-canonical partition sum reads: 

Z Mc y{ P V, V) = J W v (E n , E c , pV, V)e-^ E ^ Ec dE a dE c , (6) 

where Wv{E Q , Eq, pV,V) is the density of states. For any finite value of /3c, we treat 
the system as equilibrated at the temperature T = l//3 n and described by an equivalent 
effective charge q^ s = \(3c/(3 n . When /3 C = 0, the Coulomb energy E c does no more 
influence the partition sum and the system behaves as uncharged, reducing to the standard 
Ising model. The generalised grand potential is defined as 

Zl, Pc , a {V) = J Z Mc (pV, V)e-^ v d P , (7) 

where the fugacity is linked to the chemical potential by a = —(3 n p. When (3q = /3 n the 
ensemble coincides with the conventional grand-canonical form. 

We sampled the density distribution from eq.(J7j) for a chemical potential value p = 
p c = 3e, by employing a standard Metropolis technique [E]. An example is presented 
in Fig. [U where, for a fixed value of /3c, several density distributions are calculated for 
different temperatures T = l//3 a . All bimodal distributions present a peak for the gas 
phase (low-density) and one for the liquid (high-density) phase. Their height is equal 
because of the choice of p = p c . All distributions are composed together to draw a phase 
diagram: the ridge of the overall distribution of density and temperature delimits with 
its gas branch and liquid branch the phase coexistence. The two branches join at the 
limiting temperature. 

4. Results 

For two cases, either (3 n = (corresponding to the standard Ising model), or /3 n = Pc, 
and for a series of different lattice sizes L, we collected a series of calculations of the phase 
diagram. The results are schematized in Fig. |21 The phase-coexistence region in the 
system (3 n = shrinks for progressively larger lattice sizes according to the characteristic 
finite-size scaling behaviour of the Ising system. In particular, since the asymptotic value 
of the limiting temperature Ti im = lim^oo Ti im (L) corresponds to a critical point, finite- 
size scaling requires that Tn m (L) evolves as T[i m (L) — T\ im oc L~ x l v [US]. This is comfirmed 
by the calculated values of the critical exponent v and the critical point Xn m , which 
correspond to the Ising values. 

At variance with the Ising behaviour, the phase diagram of the system /3 n = (3q does 
however not vary with the lattice size. This first sign of incompatibility of the frustrated 
system with the Ising model, is also a first indication of quenching of the critical behaviour. 
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v rules in fact the divergence of the correlation length at the critical point, which evolves 
as £ oc t~ u with t = T/T lim — 1. If a critical point existed, it would correspond to a very 
large value for v, which would then be hardly compatible with a diverging correlation 
length. 

The major result we infer from comparing the two systems and the corresponding 
asymptotic limiting temperature, is that the coexistence region expands under the action 
of a Coulomb field. 

The connection between the increase of the limiting point and charge fluctuations can 
be searched in the evolution of the event distribution with temperature and Coulomb-field 
strength. Fig. illustrates the evolution of the limiting temperature with the strength 
of the Coulomb field and the phase coexistence is indicated for the uncharged system 
Pc — and the frustrated system fi c = /3 n , in correspondence with the phase diagrams 
shown in Fig. |21 We analyzed the event distribution along the path 1-2-3-4, determined 
by four different thermodynamic situations. 

In Fig. the distributions are shown, as a function of the two Hamiltonian compo- 
nents, (1/2V) Ydj' n i n j an d CV2V) SLj 9i?j/ r ij- Point (1) belongs to the region of phase 
coexistence sampled for the uncharged system fi c = and presents a bimodal pattern 
with respect to the abscissa, which corresponds to the nuclear-energy density E n /V, so 
that pure liquid or gas partitions are preferred to mixed events. Point (2) corresponds 
to the limiting temperature of the uncharged system /3c = 0, where mixed partitions are 
favoured. To pass from point (2) to point (3) the temperature is left unchanged while flc is 
increased from till it equals /3 n . Once introduced the Coulomb field, (A/2V) Ya^- s <?i<?j/ r 'ij 
represents the Coulomb-energy density Eq/V, which should be minimized in the frus- 
trated system. In order to satisfy this constraint, the positive-charge distribution searches 




Figure 2. Left panel. Calculations of phase diagrams for the systems (3 n = (Ising- 
like) and (3 n = j3c (frustrated). Calculations for a series of different lattice sizes L are 
associated to different limiting temperatures for the system (3 n = 0; they overlap for the 
system (3 n = [3c- Right panel. For the system f3 n = the test of the finite-size scaling 
law Tli m (L) — T u m oc L~ x l v gives Ising values for the critical exponent v and the critical 
temperature T[ im . No scaling is observed for the system /3 n = A> this is a first indication 
of the quench of criticality in presence of a Coulomb field. 
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Figure 3. (A) Increase of the limiting point and expansion of the phase-coexistence region 
with increasing frustration (i.e. for higher (3q/ (3 n ratios). Four points determine different 
thermodynamical situations: points (1) and (3) belong to the coexistence region, points (2) 
and (4) are at the limiting temperature; while points (1) and (2) refer to the system f3c = 
(Ising-like), points (3) and (4) refer to the system @ c = (3 n (frustrated). The calculation 
refers to L — 10; (B) In correspondence with the points marked in panel (A), the prob- 
ability distribution is evaluated for the two Hamiltonian components, (1/2V) Sy^i^jj 
which corresponds to the nuclear-energy density E Q /V, and (A/2V) 9iQj/ r ij ; repre- 
senting the Coulomb-energy density E c /V in the case of the frustrated system. Contour 
plots refer to L = 10; logarithmic cluster plots refer to L — 20. Points (1) and (3) manifest 
bimodal patterns. The critical scaling effect exhibited at point (2) quenches at point (4). 
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uniform patterns, which are better compensated by the uniform background of negative 
charges. As a result, the system moves to the phase-coexistence region to search pure 
partitions, and bimodality appears. In the frustrated system, the temperature should be 
further increased with respect to the limiting temperature of the system (3q = 0, point 
(3), in order to find the limiting temperature, as indicated at point (4). 

Fig. |33 presents the event distributions sampled for two different lattice sizes, in order to 
test scaling effects which could indicate the presence of critical points. At a critical point 
the correlation length £ diverges and determines the exponential decay of the correlation 
function c(rij) = (^i^j) — ( n i)( n j) according to the expression er(r) oc e~ r ^ ■ r~^ D ~ 2+r, \ 
where D is the space dimension and rj is a critical exponent. On the basis of the properties 
of the Cij matrix, cr(r) is related to the mean Coulomb-energy density by 



It is therefore typical of a critical point to exhibit a divergence of the quantity <?i?j Ay) 
for progressively larger lattice sizes. This indicates that, at point (2), the limiting point 
is compatible with a critical point, as we already expected from the finite-size scaling law 
tested in Fig. |2] for the critical exponent v. On the contrary, no scaling is observed at the 
point (4). Indeed, as indicated by Eq.(jHJ), when the Coulomb field is introduced in the 
system, the diverging quantity would be the Coulomb-energy density E c /V. To avoid 
such a singularity, the correlation length keeps a finite value at the limiting temperature 
and the critical character of the limiting temperature is lost A more formal discus- 

sion on the suppression of criticality on the basis of finite-size-scaling analysis is presented 
in ref . [ El ■ 

5. Conclusions 

By employing a specific model for studying the thermodynamical features of a frustrated 
system, we could define general properties of a neutral system in presence of charge 
fluctuations, in analogy to compact-star matter. The study of the phase diagram and the 
phase-transition phenomenology indicated that the introduction of a Coulomb field has 
the effect of increasing the limiting temperature without preserving any critical character. 

This result is the opposite as found for many other physical systems subject to Coulomb 
frustration in the absence of global charge neutrality, from hot atomic nuclei to frustrated 
ferromagnets [IHinBlEIIlEIllII]- This discrepancy with the phase-transition phenomenology 
characteristic of atomic nuclei, indicates that a connection between compact-star matter 
and nuclei is not trivial. Conversely a widening of the density range connected to the 
pasta phases was found [ ITT] within the RMF model, when the Coulomb field is included 
under the constraint of global charge neutrality over the Wigner-Seitz cells. Such a finding 
is directly compatible with the thermodynamic phenomenology we are drawing. 

It was also discussed [ ELI that neutrino propagation can be highly concerned by 
the presence of pasta phases and characterizes the process of neutron-star formation and 
cooling. With respect to homogeneous phases, the inhomogeneous phases (like pasta 
phases) present in the coexistence region are characterized by a larger opalescence to 
neutrino propagation. In the presence of a critical point the opalescence would diverge as 




(8) 
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a direct consequence of the divergence of the correlation length, and neutrinos would be 
trapped. On the contrary, we conclude that the medium stays "grey" at any temperature. 
This result seems also connected to the observation of small opacity for the transport of 
neutrinos discussed in ref . [ 122] • 
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